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Abstract 

Differential cross sections for electron collisions with the O2 molecule in its ground X^J,~ state, as well 
as excited a^Ag and b^'^g states are calculated. As previously, the fixed-bond R-matrix method based on 
state-averaged complete active space SCF orbitals is employed. In additions to elastic scattering of electron 
with the O2 X^'^g, a^Ag and b^'^g states, electron impact excitation from the X^l,g state to the a^Ag and 
b^'^g states as well as '6 eV states' of c^l,~, A'^A„ and A^l,'^ states is studied. Differential cross sections for 
excitation to the '6 eV states' have not been calculated previously. Electron impact excitation to the b^'^g 
state from the metastable a^Ag state is also studied. For electron impact excitation from the O2 X^'^' state 
to the b^'^g state, our results agree better with the experimental measurements than previous theoretical 
calculations. Our cross sections show angular behaviour similar to the experimental ones for transitions 
from the X^l.^ state to the '6 eV states', although the calculated cross sections are up to a factor two larger 
at large scattering angles. For the excitation from the a^Ag state to the b^'^g state, our results marginally 
agree with the experimental data except for the forward scattering direction. 
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I. INTRODUCTION 



A detailed knowledge of electron collisions with the oxygen molecule is important for the 
physics and chemistry of both laboratory and astrophysical plasmas. In particular, recent attempts 
to understand the electrical discharge oxygen-iodine laser have suggested that excited electronic 
states of O2 molecule play an important role lliyj]. In a previous paper Is^ (henceforth denoted I), 
we studied integral cross sections for electron collisions with the O2 molecule in its excited a^Ag 
and b^llg states, in addition to the much studied electron scattering by the O2 ground state. 
We used the fixed-bond R- matrix method with 13 target states represented by valence configuration 
interaction wave functions. State-averaged complete active space SCF (SA-CASSCF) orbitals 
with Gaussian type basis functions were employed. The calculated cross sections for electron 
impact excitation from the a^A„ state to the b^l,'!: state at 4.5 eV agree well with the available 

n 

experimental data of Hall and Trajmar ||4|]. Although elastic scattering of electrons by the a^Ag 
and b^l,g states was also studied, we could not find any experimental data for comparison. 

In I, theoretical and experimental integral cross sections were compared. However, differential 
cross sections (DCSs) provide a more stringent test of theory and are often easier to measure 
reliably than integral cross sections. For electron impact electronic excitations, calculations which 
give good integral cross sections often give DCS's which differ significantly from those observed 
experimentally. In this paper, we present DCSs for the corresponding processes calculated using 
the same R-matrix model. 

Previous experimental and theoretical studies in the field of electron O2 collisions are well 
summarized by Brunger and Buckman Is]. Here we only review works relevant to this paper. The 
DCSs of electron collisions with O2 molecule have been measured by many experimental groups. 
In particular, electron impact excitations to the low-lying a'Ap and Z?'Sl_states have been studied 
experimentally by Trajmar et al. 10], Shyn and Sweeney fl, Allan Q], Middleton et al. fl, 
and Linert et al. [10]. In contrast to these experimental works, only Middleton et al. lO] report 
calculations of DCSs for these excitation processes. Some of the more recent measurements have 
focused on electron impact excitations from the X^'^~ ground state to the '6 eV states', i.e., the 
c^Z^, A'^Aj, and A^S^ states which are also called the Herzberg pseudocontinuum llli lisl . 
Although these DCSs are not state-resolved in most case, Shyn and Sweeney IibI obtained cross 
sections for excitation to the individual electronic state within the '6 eV states'. In this paper, we 
also calculate the DCSs of this process using the fixed-bond R-matrix method, since no previous 
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theoretical calculation exists. Up to now, there is only one measurement of DCSs for electron 

n 

collisions with electronically excited O2 molecule. Hall and Trajmar |4] obtained excitation cross 
sections from the O2 to the b^llg state at electron impact energy of 4.5 eV. Their integral cross 
section was compared with our R-matrix results in I. 

In this paper, details of the calculations are presented in section 2, and we discuss the results 
in section 3 comparing our results with previous theoretical and available experiments. Then the 
summary is given in section 4. 



II. THEORETICAL METHODS 

The R-matrix method itself has been described extensively in the literature ll^fislfl^ as well 
as in I. Thus we do not repeat a general explanation of the method here. We used a modified version 
of the polyatomic programs in the UK molecular R-matrix codes [ 14] to extract T-matrix elements 
of the electron O2 scatterings. These programs utilize Gaussian type orbitals (GTO) to represent 
target molecule as well as a scattering electron. Although most of the past R-matrix works on 
electron O2 collisions have employed Slater type orbitals (STO), we select GTOs mainly because 
of the simplicity of the input and availability of basis functions. The SA-CASSCF orbitals are 
imported from the calculations with MOLPRO suites of programs jnl . The use of SA-CASSCF 
orbitals improves the vertical excitation energies of the O2 target states compared to the energies 
from HF orbitals as shown in I. These target orbitals are constructed from the the [5s,3p] con- 
tracted basis of Dunning iisll augmented by a d function with exponent 1.8846, as in Sarpal et 
al. Q. Li the R-matrix calculations, we included 13 target states; X H ,a Ag, S^,c 2„,A' A„, 
A^'L-^,,B%^,VA„ l^Ug,l^Ug,VU^ and Pn„, where the last 4 n states were not included 

in the previous R-matrix studies performed by other groups. These target states were represented 
by valence configuration interaction wave functions constructed from the SA-CASSCF orbitals. 
In our fixed-bond R-matrix calculations, these target states are evaluated at the equilibrium bond 
length 7? = 2.3 ao of the O2 ground electronic state. 

The radius of the R-matrix sphere was chosen to be 10 ao in our calculations. In order to rep- 
resent the scattering electron, we included diffuse Gaussian functions up to / = 5 with 9 functions 
for I = 0,1 functions for / = 1 - 3 and 6 functions for / = 4 and 5. The exponents of these diff"use 

in the UK R-matrix codes. Details of the 
fitting procedure are the same as in Faure et al. ll2flll . In addition to these continuum orbitals, 
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we included 8 extra virtual orbitals, one for each symmetry. The construction of the 17 electrons 
CSFs for the total system is the same as in I. The R-matrix calculations were performed for all 8 
irreducible representations of the D2/, symmetry, Ag, B2u, Bju, Big, B^,, Bi,g, B2g and A„, for both 
doublet and quartet spin multiplicity of the electron plus target system. 

The DCSs are evaluated from the T- matrix elements obtained by the R-matrix calculations. As 
in Gianturco and Jain [21] and Malegat ll22ll , the DCS is expanded using the Legendre polynomials, 



(1) 



where i and j denote the initial and final electronic states of the target, respectively. In exactly the 
same way as in Malegat but for symmetry instead of Doo/, symmetry in her paper, we can 
derive an expression of the the expansion coefficients A^., which is 
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Details of the derivation are given in the Appendix. In equation El 



etc. are 3 j 



h /' k 
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coefficients, ki is the wave number of the incident electron, 5, is the spin quantum number of the 
initial target state, while S is the spin quantum number of the total system and Ms is the projection 
of the total spin. The indices T and T' run over the 8 irreducible representations of the Dih point 
group, since we employ the polyatomic version of the UK R-matrix code. The angular quantum 
numbers of the scattering electron, Z, and m, etc. in the T-matrix element T^ji^^,,^. specify the 
real spherical harmonics 5"' instead of complex form 7"', because the 5"' transform as irreducible 
representations under symmetry. This means, there are transformation matrix elements C^_,„_. 
etc. in the expression for in order to convert the index of the scattering electron from the 5'" 
representation to the 7"' representation. An expression of the matrix elements CA,m is given in 
the Appendix. Finally, we note that the summations over (Y,li,mi) should satisfy the symmetry 
relation, Y = IR(z) x IR^S^' j, with IR(z) and IR^S^'J each being an irreducible representation of 



the D2h group corresponding to the /th target state and the real spherical harmonic 5)"' , respectively. 
This relation also holds for (Y,lj,mj), {Y' ,r^,m'^) and (r',rj,m'j). 

In equation 121 the T- matrix elements belonging to different total symmetries are multiplied to- 
gether. Since these matrix elements come from different calculations, overall phases of molecular 
orbitals and target CI vectors underlying these matrix elements may be inconsistent (see 923 ). 
which may result in erroneous relative signs of these T-matrix elements. To avoid this inconsis- 
tency, we saved reference target CI vectors from the first calculation, Ag symmetry for example, 
and then aligned the overall phases of the target CI vectors in other calculations, B2u, B^u, B\g, 
Bi,g, B2g, Au symmetries, according to this reference. The same set of molecular orbitals was used 
in all these calculations. 



ni. RESULTS AND DISCUSSION 

A. Electron collisions with the O2 ^^~g ground state 

Figure [T] shows DCSs for elastic electron scattering from the O2 ^^^g state compared with 
previous theoretical and experimental results. Our results are very similar to the previous R-matrix 



cross sections of Woste et al. 112411 . The cross sections of Machado et al. 112511 were calculated 
using the Schwinger variational iterative method combined with the distorted- wave approximation. 
Their results at 5 eV are much lower than the R-matrix results at low scattering angle below 50 
degrees. Our results agree reasonably well with the experimental cross sections at 10 eV, including 



the recent results of Linert et al. 112611 for backward scattering. At 5 eV, our model significantly 
overestimates the cross sections for forward scattering compared to the experimental values. For 
example, our result is twice as large as the experimental values at 10°. This situation is the same in 
the previous R-matrix calculation of Woste et al. [24]. As discussed by Woste et al., this deviation 
can be attributed to a lack of long-range polarizability in the scattering model. For example, 
Gillan et al. [27] introduced polarized pseudostates to account for the long-range polarizability 
in electron-N2 scattering and reduced the cross sections by 50% in the threshold energy region. 



The interaction potential of Machado et al. ll25ll includes the correlation-polarization term based 



on free-electron-gas model. Probably the polarization introduced by this term is responsible for 
their better agreement with experiment at 5 eV. Since we are interested in electron collisions with 
the excited electronic states of O2 in this work, we chose not to pursue precise accuracy further 
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for the ground state elastic scattering. However, we have to be mindful that similar long-range 
polarizability problems may exist in the other low-energy electron scattering processes, especially 
elastic electron scattering of the a^Ag and b^l,g state O2 molecules, which will be discussed below. 

The DCSs for excitation to the a^A^ state at electron impact energy 5 and 10 eV are compared in 
figureElwith the previous theoretical calculation and the experimental measurements of Middleton 
et al. |19|], Shyn and Sweeney ||7|], Allan jsl and Linert et al. Iioi . The cross sections at 5 eV agrees 
well with the previous calculation and experimental data below 120°. However, our results are 
much lower than the previous calculation of Middleton et al. at scattering angle above 130°. At an 
electron scattering energy of 10 eV, our cross section deviates further from the previous calculation 
of Middleton et al. 1^] especially at scattering angle below 60° and above 140°. In contrast to the 
backward enhanced cross sections of Middleton et al., our DCSs have a slightly forward enhanced 
character. Our results agree better with the experimental data at low scattering angles than the 
previous calculation of Middleton et al. ||9 K At large scattering angle, our results deviate from the 
experimental data of Shyn and Sweeney IlZD, but agrees rather well with the recent measurement 
of Linert et al. Iioll . The R- matrix model of Middleton et al. [9] included the lowest 9 O2 target 
states and / = - 5 scattering electron orbitals with cr, n and 6 symmetry. In this work, we included 
13 target states and all components of / = - 5 scattering electron orbitals. In addition to these 
differences, Middleton et al. used HF/STO orbitals where we employed CASSCF/GTO orbitals. 
We carried out a test calculation with / = - 3 scattering electron orbitals and got almost the same 
cross sections as in Z = - 5 case, which suggests that difference in the number of target states may 
be important for the shape of these excitation cross sections. 

Figure El compares DCSs for excitation to the b^l,g state at electron impact energy of 5 and 10 
eV with the previous R-matrix calculation and the experimental measurements of Middleton et al. 
lO], Shyn and Sweeney |Q] and Allan [8]. Transitions between S+ and 2~ target states are forbidden 
at scattering angles of 0° and 180°, because the scattered electron wave function vanishes in the 
plan e defined by incident electron beam and the molecular axis for any orientation of the molecule 
J28[ ligj . As a consequence, the DCSs decrease to be zero toward and 180°. As is apparent 
from figureEl our cross sections become zero at and 180°, which is consistent with this selection 
rule. Compared to the previous R-matrix calculations of Middleton et al. [O], our cross sections 
have similar profile, but with slightly smaller magnitude at all scattering angles. Agreement with 
experiment is good at 5 eV below 120°, although our results underestimate the experimental cross 
sections at larger scattering angles. At 10 eV, the magnitude of the experimental cross sections of 
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Middleton et al. [Q] and Shyn and Sweeney |Q] do not agree with each other, however our cross 
sections are closer to the results of Shyn and Sweeney at low scattering angles below 50°. Between 
60° and 90°, our results are closer to the results of Middleton et al. 

Figure m shows DCSs for excitations to the '6 eV states' for electron impact energies of 10 and 
15 eV. Here the '6 eV states' means the group of the O2 c'S^^, A'^A„ and A^2^ states. The cross 
sections shown in figure HI are a sum of individual excitation cross sections of these 3 electronic 
states, in line with most experimental measurements. The figure includes the recent experimental 
cross sections of Campbell et al. [ H^, Green et al. jvl] and Shyn and Sweeney lisl . The individual 
cross sections are shown in figure |5l and |6l for impact energies of 10 and 15 eV, together with the 
state-resolved experimental cross sections of Shyn and Sweeney lll3ll . Our summed cross sections 
given in figure HI are backward-enhanced for both the 10 and 15 eV cases, in accordance with 
the experimental cross sections Campbell et al. [ llil and Shyn and Sweeney IibI . However, the 



forward enhancement of the DCSs at 10 eV observed by Green et al. 111211 is not reproduced by 
our calculation. The individual cross sections in figure |5] and |6l show similar angular behaviour 
compared to the experimental results of Shyn and Sweeney. However our DCSs for excitation to 
the A'^A„ state is more steep toward backward direction. Also the peak in the A^S^^ state cross 
sections is more pronounced in our calculation. Note that our results for the A^Z^ state become 
zero at and 180° as dictated by the 2" - selection rule. 



B. Electron collisions with the O2 a^Ag and b^'Eg excited states 

The DCSs for elastic electron scattering with the excited O2 a^Ag and b^llg states are shown in 
figure |71 We cannot compare them with previous theoretical or experimental work, since there is 
no available data. These DCSs show strong similarity with those of the elastic electron scattering 
with the X^'L'^ ground state in figure [H The magnitude of these cross sections is almost the same 
for the 10 eV case. All of them have a large forward peak at 0°, a small rise in the cross sections 
at 180°. The location of the minimum moves inward from 140° to 90° as the electron scattering 
energy increases. This similarity is also reflected in the integral cross sections for elastic electron 
collisions with the X^l,~, a^Ag and b^l,^ states. The profiles and magnitudes of the integral cross 
sections are basically the same for all these 3 electronic states as shown in I. The main configu- 
ration of these 3 electronic states has the form (core);r^;7r^, and this may be responsible for this 
similarity. Our R-matrix calculations tend to overestimate the elastic scattering cross sections of 
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the X^'L' state at low scattering angles, below 50°, compared to the experimental data. Consider- 
ing the strong similarity of the cross section profiles for elastic scattering from excited states and 
the ground state, our calculations may also overestimate the cross section at low scattering angle 
at low electron impact energy. 

In table U we show momentum transfer cross sections for electron elastic scattering by the 
X^^g, a^Ag and b^l,g states. As a consequence of the similarity in DCSs, the momentum transfer 
cross sections have a similar magnitude. Compared to the experimental data of Shyn and Sharp 
1 3011 and Sullivan et al. [31], our calculation overestimates the X^'Z' state momentum transfer 
cross section at 2 eV by 20 %, but underestimates the cross section at 10 eV by 8 % of Sullivan 
et al.'s value or 29 % of Shyn and Sharp's value. Our momentum transfer cross sections for the 
a^Ag and states may similarly be overestimates or underestimates depending on the electron 
impact energy. 

Figure [5] shows DCSs for electron impact excitation from the a^Ag state to the b^'Zg state. The 

figure also includes the experimental data of Hall and Trajmar lO] at impact energy of 4.5 eV. 

Our cross section profiles have characteristic features of minima around 10° and 90° and maxima 

Jl 

around 50° and 150°. They agree with the experimental cross sections of Hall and Trajmar ||4(] 
within their error bars except at 20° and 30°. The cross sections of Hall and Trajmar appear to 
increase from 50° to 0° whereas our cross sections decrease from 60° toward 10°. In the 50° - 
140° angular region. Hall and Trajmar's cross sections vary less than ours. However, a precise 
comparison is difficult because of large error bars and lack of other experimental data. 



IV. SUMMARY 



We have calculated differential cross sections for electron collisions with O2 molecule in its 
ground X^l,~ state, as well as excited a^Ag and b^llg states. As in our previous work, we employed 
the fixed-bond R- matrix method based on state-averaged complete active space SCF orbitals. In 
additions to elastic scattering of electron with the O2 X^l,~, a^Ag and b^l,g states, we studied 
electron impact excitations from the X^H' state to the a^A^ and b^llg states as well as '6 eV states' 
of c'S" A'^A,, and A^Z;*' states. DCSs for the excitations to the '6 eV states' were not calculated 
previously. We also studied electron impact excitation to the state from the metastable a^Ag 
state. For electron impact excitation from the O2 X^l,g state to the b^llg state, our results agree 
better with the experimental measurements than the previous theoretical cross sections. Our cross 
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sections show similar angular behaviour the to experimental ones for transitions from the X^l.' 
state to the '6 eV states'. For the excitation from the a^Ag state to the b^'Z^ state, our results 
marginally agree with experimental data except for the forward scattering direction. 
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FIG. 1: Differential cross sections for elastic electron collisions with the O2 X^^g state. Panel (a): electron 
impact energy of 5 eV and (b):10 eV. Thick full line represents our result. For comparison, previous theo- 
retical results of Woste et al. f24\ and Machado et a l. Ii25i1 are shown as thin lines. Symbols with error bars 
indicate experimental cross sections of Woste et al. 1241] . Sulhvan et al. |31], Green et al. |32] and Linert et 
al. 
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FIG. 2: Differential cross sections for electron impact excitation from the O2 X^'^g state to the a^Ag state. 
Panel (a): electron impact energy of 5 eV and (b):10 eV. Full line represents our result. For comparison, we 
include previous theoretical and experimental cross sections of Middleton et al. |9], experimental results of 
Shyn and Sweeney |7|], Allan and Linert et al. llflfl. 
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FIG. 3: Differential cross sections for electron impact excitation from the O2 X "Zg state to the b^'^g state. 
Other details are the same as in figure |2l 
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FIG. 4: Differential cross sections for excitation from the O2 X^'^g state to the '6 eV states' which consist 
of the c^Ej^, A'^Ak and A^E^^ states. The cross sections shown here are sum of the individual cross sections 
of these 3 states. Panel (a): electron impact energy of 10 eV and (b):15 eV. For cornparison, we include 
experimental results of Green et al. Il2ll . Campbell et al. Illi and Shyn and Sweeney lisi . 
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FIG. 5: Differential cross sections for excitation from the O2 X^E state to the individual state of the '6 



eV states'. Panel (a) shows excitation cross sections for the c'S^^ state, panel (b) is for the A'^A„ state and 
panel (c) is for the A^E^^ state. The electron impact energy is 10 eV. The experimental data was taken from 
Shyn and Sweeney lisll . 
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FIG. 6: The same as in figure |5l but for an electron impact energy of 15 eV. 
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FIG. 7: Differential cross sections for elastic scattering of the O2 excited states, (a): the Ag state, (b): the 
fe^Xp state. Each line corresponds to a different electron impact energy as shown in the legend. 
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FIG. 8: Differential cross sections for electron impact excitation from the O2 tJ^A^ state to the b^'^g state. 
Each line corresponds to a different electron impact energy, as shown in the legend. Experimental cross 
sections at 4.5 eV of Hall and Trajmar 0] are also shown. 
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TABLE I: Elastic momentum transfer cross sections in unit of 10 '^cm^. 



Electron impact enerj 


ry(eV) 


2.0 


4.0 


6.0 


8.0 


10.0 


14.0 


This work 


X% 


8.18 


6.80 


6.73 


6.06 


6.04 


6.96 




a^Ag 


8.16 


6.69 


6.57 


5.90 


5.86 


6.79 






7.95 


6.57 


6.29 


5.83 


5.85 


6.73 


Shyn and Sharp [301 




6.7 








8.4 




Sullivan et al. [31] 




6.5 


6.0 




6.2 


6.4 





20 



APPENDIX: DERIVATION OF EQUATIONH 



Derivation of the DCS formula in equation |21 is similar to that of Malegat ll22ll except for the 
use of real spherical harmonics S"' employed in the polyatomic version of UK molecular R-matrix 
codes instead of complex form Y'". For convenience of reader, brief derivation of the formula is 
given in this appendix. In the expressions below, we follow the notation of Malegat I22I1 . 

The scattering wave function describing collision of an electron plane wave with a molecule is 
expressed as, 

{x[, .., x'j,, a', r) = {x[, .., x'^)xi^^^^ {cr') e''-' + J] T,- {x[, .., (cr') Fjj (r) e''^'/r. 

(A.l) 

Here x denotes the space and spin coordinates of the molecular electrons. The primed coordinates 
refer to the molecular frame with z'-axis along the molecular symmetry axis, and the unprimed 
coordinates to the laboratory frame with the z-axis along the incident electron beam. The incident 
electron has wavenumber ki with spin projection m^, . The index i represents quantum numbers of 
the electronic state of the target molecule, F,, 5; and Ms^, whereas the index / refers to (z, m^,) 
collectively. 

In order to expand the wave function in equation lA.il a symmetry adapted N+1 -electron wave 
function is prepared as 

Tg* (x'„ ... a; r) = T-' «, ... a')S':; (f /-■. (^.2) 

where F, S and M stand for symmetry of the N+l-electron system, i.e., an irreducible represen- 
tation of the D2h group in this work, spin quantum number and its projection to the symmetry 
axis. The orbital angular momentum of the scattering electron and its projection are represented 
by /, and m,. In case of m > 0, the real spherical harmonics 5"' is related to the complex form of 
spherical harmonics Y'p as 113 3h 

^[(-ir (-1)'"/ 
V2[ 1 -i 

In the m = case, we only have 7*' and 5*' and thus the matrix element is 1. Note that 5"' behaves 
as an irreducible representation under D2/, symmetry operations, whereas Y'" does not. The spin 
coupled function in equation lA.2l is given by 

^5^^- {x\, .., x'^, cr') = J] {SjMs, \m,^\SMs)^j {x\, .., x'^)xim^^ , (A.4) 



y-ni 



(A.3) 
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where {S jMs j,\ms-\S Ms) refers to the Clebsch-Gordan coefficients and / represents F,- and 5/ 
The radial functions in equation |a21 are related to the S -matrix in the asymptotic region by 

1 



lim (r) = — 



ilimijljntj 



(A.5) 



Expanding equation lA.ll in the symmetry adapted functions of equation lA.2l gives 



7/,m, rsMi- 

By comparing the ingoing parts on the right and the left hand side, we obtain the expansion coef- 
ficient 

a™ = ^^^3 ^ Z^S ^ ^""^y^ C,,^.(S,Ms., \m,\SMs), (A.7) 

where Hy^^ ^, (a/Sy) is the rotation matrix with the Euler angles (a,/3, y) representing rotation of the 
laboratory frame to the molecular frame. The matrix element defined in equation [A.3l relates 
the spherical harmonics Yf and S"'. The collision amplitude can then be obtained by equating the 
outgoing parts, 

Fij(r) = 2 Z ^""^^ ^\ '--'^^'(SiMs„ ^m,,\SMs){SjMs^, ^m,^\S Ms) 

limiljnij TS Ms Afiv V ' ^' 

XD'- 1 (apy) di, ia/3y) (r) C,,„,C;,„T™^.^^. (A.8) 

Here we use the T-matrix elements T^.f^l, instead of the S-matrix. 

By summing over the final states and averaging over the initial states and the molecular orien- 
tation (q',jS, y), the differential cross section is expressed by the Legendre polynomials expansion 
[Uwith expansion coefficients given by equation |2l 
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